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q ! Abstract 

In the inflation models, the relic gravitational waves (RGW) generated in the inflation stage, 
and evolved in the Universe until now. In the different cosmological evolution models, one 
can get different gravitational waves power spectrum. In this paper, we give a simple formula 
to estimate this spectrum in a general cosmological model. From this formula, one can easily 
find the relation between this power spectrum and the cosmological evolution models. The 
spectrum includes all the information about the evolution of the scale factor a from inflation 
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to now. So RGW is a more clear fossil, which records the cosmological evolution information, 
and be a useful complementarity for another fossil— Cosmic Microwave Background Radiation 
(CMB). 
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In the inflation expansion cosmic models, the stochastic background of relic gravitational waves can 
generate from the quantum fluctuate, which is important in cosmology and has been extensively studied 
in the past [I]. The power spectrum of it, which is observed today, depends not only on the details of the 
early stage of inflationary expansion, but also on the expansion behavior of the subsequent stages [2, 3], 
especially, some stages as the pre-big bang stage, the reheating stage, the cosmological QCD transition, 
the e + e~ annihilation process, the accelerating expansion process, which we did not know very clearly in 
physics until now. Their information has also been recorded in the RGW power spectrum. Most of these 
stages occurred before the CMB photon decoupled, which makes it difficult to study them from CMB. 
So research on RGW power spectrum is nearly the only very way to study this unknown process of the 
Universe. In this paper, I will study this power spectrum in a general cosmological model, and find that 
there is a very simple formula to estimate it. From this formula, one can easily find the relation between 
the spectrum and the evolution of the Universe. Especially, the spectrum index directly relates to the 
state equation of matter which dominated the evolution of the Universe at some stage. This gives a simple 
way to estimate the cosmic evolution process from the RGW power spectrum. 



In the spatially flat Robertson- Walker spacetime, the metric which includes the gravitational waves is 

ds 2 = a 2 {r)[dT 2 - (5ij + hjjUiyd-r 1 ]. (1) 

where r is the conformal time, the perturbations of spacetime h^ is a 3 x 3 symmetric matrix. The 
gravitational wave field is the tensorial portion of hij, which is transverse-traceless dih^ = 0,5^ = 0. 
Since the relic gravitational waves are very weak, \hij\ <C 1, so one needs just study the linearized field 
equation: 

0„(>/=^Mx,t)) = O. (2) 

In quantum theory of gravitational waves, the field hij is a field operator, which is written as a sum of 
the plane wave Fourier modes 

M*,r) = ^fi± r ^(kJ^IaWfcWW^ + of^fWe-^, (3) 
(27r)2 t^i J -°° v2fc 

where Ipi = \[G is the Planck's length, the two polarizations e^(k), A = 1,2, are symmetric and 
transverse-traceless /c*e^(k) = 0, <5^e^(k) = 0, and satisfy the conditions e^ lJ (k)e^(k) = 25\y, 
and €-ij\— k) = e-^(k), the creation and annihilation operators satisfy [a^\a\^, X ■*] = 5\\>5 3 (k — k'). 

For a fixed wave number k and a fixed polarization state A, the wave equation reduces to the second- 
order ordinary differential equation [3, 4] 

4 A) "(r) + 2^«'(r) + ^„W=0, (4) 

where the prime denotes d/dr, and k = |k|. Since the equation of h^\r) for each polarization is the same, 
we denote ti£\r) by /ifc(r) in below. 

The power spectrum h(k, r) of relic gravitational waves is always defined as 

h(k,r) = ^Prk\h k (r)\. (5) 
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The spectral density parameter Q g (v) of gravitational waves is defined through the relation p g /p c = 
J £l g (y)^j-, where p g is the energy density of the gravitational waves and p c is the critical energy density. 
One reads 

n g (v) = ^h 2 (k,T H )(— ) 2 , (6) 

3 u H 

where th is the conformal time of now and i/jj is the Hubble frequency of now. Which can not exceed 
1CT 6 as required by the rate of the primordial nucleogenesis[5]. 

In the inflation models, RGW generated from the initial quantum fluctuate during the inflationary 
stage. In general, in the inflation stage, the scale factor has a form 

o(r) = 1 \t\ 1+i3(t \ -cxxr<0, (7) 

Normally, (3 is a constant. If (3 = —2, the inflation is a simple de Sitter expansion. But recently 
from observed results of the primordial scalar power spectrum from WMAP[6], one finds the inflation is 



a very complicated process, here we assume (3 is a slowly evolving function with time r. For a given wave 
number k, it crossed over the horizon at a time Tj, i.e. when the wave length Aj = 27ra(ri)/k is equal to 
\/H(ji), the Hubble radius at time Tj. Eq. (7) yields 1/H(ti) = £o|7"i| 2 /|1 + (3\, and for the exact de 
Sitter expansion with (3 = —2 one has Hfa) = Iq. Note that a different k corresponds to a different time 
Tj. Now choose the initial condition of the mode function hk(r) as [3] 

IM^OI = T7~T- (8) 
Then the initial amplitude of the power spectrum is 

h{k,T i ) = &^l pl H{T i ), (9) 

From Aj = 1/H(ji) it follows that = Jp. So the initial amplitude of the power spectrum is 

h(k, n ) = Ak 2+ ^ = Ak 2 ^- 1 ^, (10) 

where the constant 

A = 8V^l pl \l + (3(r i )\/l . (11) 

For the case of (3 = —2, the initial spectrum is independent of k, this is the so called scale independent 
initial power spectrum. Its amplitude directly relate the inflation Hubble constant. 

The evolution of RGW in the following stage satisfies the two order differential equation (4). Here 
we only suppose that the Universe is always a decelerating expand, (from the research of SNIa, WMAP 
and LSS, one finds that now the Universe is accelerating expansion, we have researched its effect on the 
spectrum and found it was little [3]. ) If the matter state equation satisfies p = oop, where p and p are the 
pressure and energy density of the matter respectively, (uj is a constant, for radiation, to = 1/3; and for 
matter uj = 0), one can easily find that the scale factor has a form: 

a = b[T-T 1 } a , (12) 

where a = , b and t\ are constant. For the component of the Universe vary with time, we here think 
the uj is varying slowly with time, this is a better way to describe the evolution of the Universe than the 
method which separated the Universe into several different evolution stages. So the scale factor can be 
described in the following form: 

a = b(r)[r - n(r)] a « , (13) 

where 6(r), n(r) and a(r) are all slowly evolution functions with time. This way can be used to study all 
kind of cosmic evolution models. Especially, it is useful for studying the effect of the cosmological phase 
transition on the RGW power spectrum. 

In this model, when the wavelength of RGW is much larger than the horizon, the amplitude of the 
spectrum keeps constant. Once the gravitational waves enters the horizon, the spectrum will have a 



damping along with the expansion the Universe. We assume that at a certain short stage, the scale factor 
has a form of 

a = b'[r-r[] a ', (14) 

For b'(r), r{(r) and a'(r) vary slowly with time, they can be assumed as constant in this stage. We study 
the gravitational waves which enter the horizon in this stage. A wave with wavenumber k, its amplitude 
satisfies [3] 

h(r H , k) = h(n, k)a*{T)/a(T H ), (15) 

where a*(r) is the scale factor when the wave k just enters horizon, and a(r#) is the now scale factor. 
From the scale factor formula and the relation k(r — r{) ~ 1 when wave enters horizon one find 

a*(r)/a(T H ) = Kl ^ )k \^ (16) 

where bn = b{Tn)- So we can get the spectrum 

h(k,r H ) = (A/bH^l/k)^^ 1 ^ 2 '^ = ^l P iH{r^/k)/b H )k^ TH) k-^l k \ (17) 

For the slowly varying with time of 6(t), n(r) and a(r), we can use this formula to estimate the RGW 
spectrum for all k. From this, one can also get the spectrum index 

n{k) = d\nh(k,T H )/dk ~ 0(-l/k) + 2-a(l/fc), (18) 

So the matter state equation when the horizon of the Universe equate to 1/k is 

"W ="('/*) = 3 ^ ( . 1A) 2 +2 - BW1 -'. < 19 > 

In general, j3 + 2 is nearly zero for inflation, so this formula directly relates the RGW spectrum index and 
matter state equation in the Universe. For the waves with wavelength long enough to have not reenter the 
horizon until now, they have the power spectrum of their initial condition 

h{k,r H ) = h(k,n) = 8^lpiH(n) = Ak 2+I3 ^ k \ (20) 

and 

n(fc)~/9(-l/fc)+2. (21) 

These are the approximate formulae, which can evaluate the RGW power spectrum (If one want to know 
the exact spectrum, it must depend on the numerical calculation or the complicated exact calculation), 
but it is enough for knowing the physical information which included in the RGW power spectrum. From 
formula (20), one can find these spectrum keeps its initial value, which reflects the physics of the inflation. 
Its amplitude is exactly the Hubble parameter of the inflation, its change with k is also the change of 
Hubble parameter with time. So from which we can know not only the energy scale of the inflation which 
closely relates the inflation physics origin, but also the inflation evolution information, which is also very 



important for the very damping of the CMB anisotropics power spectrum in the very large scale, and the 
large running of index of the primordial scalar perturbation observed by WMAP[6, 7], some people think 
that it reflects the complicated evolution of early inflation. 

For the waves with larger frequency, the spectrum formula (17) also includes the evolution information 
of the Universe after inflation. Its amplitude is dependent on the all evolution process of Universe after 
cosmic horizon scale became larger than the wavelength. Especially the spectrum index n s , which directly 
relates to the matter state equation when the cosmic horizon scale equate to the wavelength, and the 
variety of the spectrum index with wavelength is also the variety of the matter state equation in the 
Universe with time. Actually, we can only detect the RGW power spectrum in several narrow frequency 
span, for example the laser interference explorers [8], which is sensitive for the wave with frequency of 
1(T 3 ~ 10 2 Hz, the method using the millisecond pulsars[9], which is sensitive at 10 9 ~ 10 7 Hz, the 
method by researching the anisotropy and polarization of CMB [10], which is sensitive at 10~ 18 ~ 10~ 16 Hz. 
So even if we have detected RGW at some frequency, we also can't know clearly the inflation energy scale 
unless we know all information of Universe after the inflation stage, but we can research the matter state 
equation in a special stage of the Universe by the index of the power spectrum, which directly relates the 
matter state equation uj at the time r = 1/k by formula (19). So I think the simple RGW power spectrum 
in all frequency can reflect the all evolution information of the Universe since inflation, and it can not be 
contaminated by other information as matter perturbation, the interaction between the matters. So it is a 
very clear fossil for study the cosmic evolution, and be a useful complementarity for another fossil— CMB, 
especially for studying of the Universe stage before the CMB photon decoupled. 

References 

[1] S.Sasaki, Prog.Theor.Phys. 76,1036 (1986); A.A.Starobinsky, JEPTLett. 30, 682 (1979); A. Rubakov, 
M.V.Sazhin, and A.V.Veryaskin, Phys.Lett.115B, 189 (1982); R.Fabbri and M.D. Pollock, Phys.Lett. 
125B, 445 (1983); L.F.Abbott and M. B. Wise, Nucl.Phys. B224, 541 (1984); L.F.Abbott and 
D.D.Harari, Nucl.Phys. B264, 487 (1986); V.F.Mukhanov, H.A.Feldman, and R.H.Brandenberger 
Phys.Report 215,203 (1992); D. H. Lyth, A. Riotto, Phys.Rept. 314 1 (1999); 

[2] B.Allen, Phys.Rev.D 37, 2078 (1988); V.Sahni Phys.Rev.D 42, 453 (1990); L.Grishchuk, 
Class. Quant. Grav. 14 1445 (1997); Lecture Notes Physics 562, 164 (2001), in "Gyros, Clocks, In- 
terferometers...: Testing Relativistic Gravity in Space" , Lammerzahl et al (Eds); A. Riszuelo and J-P 
Uzan, Phys.Rev.D 62, 083506 (2000); H. Tashiro, K. Chiba, and M. Sasaki, Class. Quant. Grav. 21 
1761 (2004); A.B. Henriques, Class. Quant. Grav. 21 3057 (2004); 

[3] Y.Zhang, Y.F.Yuan, W.Zhao, Y.T.Chen, Class.Quant.Grav. 22, 1383 (2005); Y.Zhang, W. Zhao, 
Y.F. Yuan and TY. Xia Chin.Phys.Lett. 20, 1871 (2005); 

[4] L.Grishchuk, Sov.Phys.JETP, Vol. 40, No.3, 409 (1974). 

[5] E.Kolb and R.Turner, The Early Universe (Addison- Wesley, Reading, MA, 1990); B.Allen, gr- 
qc/9604033; M.Maggiore, Phys.Rep. 331, 283 (2000); 

[6] C.L.Bennett, et al. Astrophys.J.Suppl. 148, 1 (2003); D.N.Spergel, et al. Astrophys.J.Suppl. 148, 175 
(2003); 

[7] G. Hinshaw, et al. Astrophys.J.Suppl. 148, 135 (2003); H.V.Peiris, et al. Astrophys.J.Suppl. 148, 213 
(2003); 



[8] http: // www.ligo.caltech.edu; http://www.geo600.uni.hannover.de; http://http://lisa.jpl. nasa.gov/; 
LIGO Scientific Collaboration: B. Abbott, et al, astro-ph/0507254; 

[9] S. Detweiler, Astrophys. J. 234, 1100 (1979); A.N. Lommen, astro-ph/0208572; 

[10] M. Kamionkowski and A. Kosowsky, Ann. Rev. Nucl. Part. Sci. 49, 77 (1999); W. Hu and S. Do- 
delson, Ann. Rev. Astron. Astrophys. 40, 171 (2002); M. White and J. D. Cohn, astro-ph/0203120; 
S. Dodelson, Modern Cosmology (Academic Press, 2003); A. Liddle and D. H. Lyth, Cosmological 
Inflation and Large-Scale Structure (Cambridge University Press, 2000); Wen Zhao and Yang Zhang, 
astro-ph/0509753, astro-ph/0510355; 



